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Abstract 

A number of subspecies of the eastern hedgehog Erinaceus concolor has been described in the southern 
Balkan peninsula. Despite extensive morphological and karyological descriptions, their taxonomic sta¬ 
tus is still unclear. In the present study we have attempted to identify this status in different populations 
from mainland and insular Greece and to clarify their relationships to subspecies from adjacent areas. 
Multivariate analysis of 21 somatic and cranial characters was performed on 191 specimens from 
Greece, and also on published data (76 from Bulgaria, 21 from Yugoslavia and 9 from Asia Minor). 

Our analysis showed that throughout the mainland of the southern Balkan peninsula, one large¬ 
sized subspecies, E. c. bolkayi, is distributed. Two insular subspecies of smaller body size are found on 
the Greek islands: E. c. nesiotes on Crete, Kyklades islands and Ionian sea islands and E. c. rhodius in 
the eastern Aegean sea islands Rhodos, Samos, and Chios. The subspecies E. c. transcaucasicus also ap¬ 
pears to occur on the islands of Lesvos and Kos. 


Introduction 

The eastern hedgehog Erinaceus concolor has an extensive distribution in eastern Europe 
and west Asia: Poland, Czech Republic, Slovakia, Austria, north-eastern Italy, Balkan 
peninsula, Turkey, southern Russia, west side of the Caspian Sea and western Siberia as 
well as in Israel, northern Iraq and northern Iran (Corbet 1988; Holz and Niethammer 
1990). In this vast area, the hedgehog displays such a degree of variation, particularly in 
colour and size, that a considerable number of subspecies has been described. A melanic 
form Erinaceus concolor concolor (or E. c. ponticus) is distributed on the eastern shore of 
the Black Sea and a pale form E. c. palidus in west Siberia (Ognev 1928; Corbet 1988). 

Using a multivariate analysis of skull measurements, Holz (1978) showed that apart 
from the western and eastern hedgehog, a third type exists in south-eastern Europe (Ro- 
mania-Bulgaria) with a smaller cranial capacity and a longer cranium. It is not certain 
whether this type belongs to the species E. concolor, to the subspecies E. c. transcaucasicus 
(Satunin, 1905), or to a third species not yet described. 

Kratochvil (1980) states that the hedgehogs inhabiting lowland steppes of the Balkan 
peninsula have a smaller mean capacity of the skull cavity and belong to the form de¬ 
scribed as E. c. drozdovskii (Martino 1933). Thus, according to Kratochvil (1980) the 
European continent is populated by two separate subspecies of Erinaceus concolor differ¬ 
ing in somatic and craniological characters i. e. roumanicus and drozdovskii. The form 
drozdovskii has been described as a “morpha” of the subspecies Erinaceus roumanicus 
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bolkayi (Martino 1930), being lighter in colour and larger in size and distributed in the 
mountaneous regions of Montenegro, Yugoslavia (Martino 1933; Petrov 1940). Some 
authors considered bolkayi and drozdovskii two separate subspecies of the eastern hedge¬ 
hog, larger in body size than the subspecies E. c. roumanicus (Markov 1957; Markov and 
Dobrijanov 1974; Malec and Storch 1963; Dulic and Muric 1967; Ondrias 1966). Ac¬ 
cording to Wettstein (1966) the hedgehogs from Asia Minor belong to E. c. transcaucasi- 
cus (Satunin, 1905). He considered it as a subspecies synonymous with drozdovskii, which 
belongs to Anatolian populations. This subspecies has also been found on Kos (de Beaux 
1929). 



Fig. 1. Distribution map of the subspecies of Erinaceus concolor in Greece. • = E. c. bolkayi, 

A. = E. c. nesiotes, M = E. c. rhodius, * = E. c. transcaucasicus. Continental Greece: a = Attica 
(ai = Koropi, a 2 = Parnitha), b = Boiotia (Thebes), e = Evoia (ei = Aliveri, e 2 = Edipsos), 1 = Fthiotida 
(Lamia), f = Fokida (fj = Lidoriki, f 2 = Tichio), g = Aetoloakarnania (Agrinio), i = Epiros (loannina), 
s = Thessalia (si = Karditsa, S 2 = Trikala, S 3 = Kalambaka, S 4 = Elassona), m = Makedonia (mi = Fiorina, 
m 2 = Alexandria, m 3 = Kilkis, m 4 = Serrae), t = Thrace (ti = Alexandroupolis, t 2 = Soufli). Peloponni- 
sos: 1 = Kastritsi (Patra), 2 = Pyrgos, 3 = Kiparissia, 4 = Astros (Kinouria), 5 = Korinthos. Crete: 

1 = Chania, 2 = Rethymno, 3 = Iraklio, 4 = Agios Nikolaos, 5 = lerapetra. 
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An extensive karyological investigation of hedgehogs in Greece has shown that three 
subspecies of the eastern hedgehog E. concolor are distributed throughout the country 
(Giagia and Ondrias 1980). The large-sized hedgehog of the Balkan peninsula was first de¬ 
scribed in the mainland as E. c. drozdovskii. However, a recent morphological study 
showed that it should be denominated E. c. bolkayi (Giagia and Ondrias 1995). An insular 
subspecies, E. c. nesiotes (Bate, 1905), is found on the island of Crete and is karyotypically 
identical to the mainland form, while a second insular subspecies E. c. rhodius (Festa, 1925) 
is found on the island of Rhodos and differs karyotypically from the other two subspecies 
(Giagia and Ondrias 1980). Some hedgehogs from the island of Vis in Yugoslavia have 
been described under the name E. c. nesiotes (Dulic and Tvrtkovic 1979). 

In the present study, we attempted to give a global view of the morphometric relation¬ 
ships among the populations of E. concolor throughout Greece and to clarify their rela¬ 
tionships with the neighbouring subspecies that have been previously described. 


Material and methods 

The present analysis was performed on 191 specimens constituting eleven populations collected from 
the mainland of Greece and its numerous islands (Fig. 1). Only adults were used to minimize the effect 



Fig. 2. Measurements taken on skulls of Erinaceus concolor (GSL = greatest skull length, 

CBL = condylobasal length, BL = basal length, ZB = zygomatic breadth, MB = mastoid breadth, 

IC = interorbital constriction, NL = nasal length, RB = rostrum breadth, PB = postorbital breadth, 
PW = palate width, UT = upper toothrow, ML = mandible length, LT = lower toothrow, HS = height of 
skull, PL = palate length, HRM = height of ramus mandibulae, RA = rostral angle). 






















































132 


Eva B. Giagia-Athanasopoulou and G. Markakis 


of allometric variation associated with growth. Twentyone measurements, four bodily (HB = head and 
body length, TL = tail length, EL = ear length, and HFL = hind foot length) and seventeen cranial mea¬ 
surements (Fig. 2), were selected for the present analysis. In addition, two indices that have been used 
for the study of the hedgehog (Wettstein 1942; 1966; Ruprecht 1972) were calculated: maxillary index 
(Mi = a - b/a - c) and length index (Li = CBL/ZB). 

Additional data from Bulgaria (Markov 1957), Yugoslavia (Dulic and Tvrtkovic 1979; Martino 
1933; Petrov 1940) and Asia Minor (Wettstein 1966) were used for comparative purposes. 

A list of all populations including the measurements used is given in table 1. Cranial measurements 
were taken as defined by Hrabe (1976) with a precision of 0.1 mm or 0.5 degree in the case of the ros¬ 
tral angle (Ruprecht 1972). 

Data were subjected to canonical (discriminant) analysis (CA) and principal components analysis 
(PCA) after their logarithmic transformation. MANOVA was used to detect differences between popu¬ 
lations (Morrison 1976; Reyment et al. 1984) and Mahalanobis distances between population centroids 
in order to construct a UPGMA phenogram (Sneath and Sokal 1973). The statistical packages 
SYSTAT (Wilkinson 1987) and NT-SYS (Rohlf 1992) were used throughout this study. 


Table 1. Populations, sample sizes and variables (in common) studied. 


No 

Description 

Code 

Sample size 

Variables used 

1 

Continental Greece 

CG 

49 

HB, TL, HFL, EL, GSL, CBL, BL, 
ZB, MB, IC, NL, RB, PB, PW, UT, 
ML, LT, HS, PL, HRM, RA 

2 

Peloponnisos 

PE 

76 


3 

Crete {E. c. nesiotes) 

CR 

15 

” 

4 

Rhodos {E. c. rhodius) 

RH 

9 

” 

5 

Lesvos 

LE 

12 

” 

6 

Kyklades islands 

KY 

4 

” 

7 

Kerkyra 

KE 

5 

” 

8 

Kephalonia 

KP 

3 

” 

9 

Kos 

KO 

5 

” 

10 

Samos 

SA 

7 

” 

11 

Chios 

CH 

6 

” 

12 

Bulgaria {E. c. drozdovskii) 

Bd 

23 

HB, TL, HFL, EL, GSL, CBL, ZB, 
MB, IC, NL, RB, PB, UT, ML, LT 

13 

Bulgaris {E. c. bolkayi) 

Bb 

23 

” 

14 

Bulgaria {E. c. roumaniciis) 

Br 

30 

” 

15 

Yugoslavia {E. c. drozdovskii) 

Yd 

6 

HB, TL, HFL, EL, CBL, ZB, UT, 
ML,LT 

16 

Yugoslavia {E. c. bolkayi) 

Yb 

6 

HB, TL, HFL, EL, CBL, ZB, MB, 
NL, RB, UT, ML, LT 

17 

Yugoslavia {E. c. nesiotes) 

Yn 

9 

” 

18 

Asia Minor (E. c. transcaucasicus) 

Tt 

9 

CBL, ZB 


Results 

Preliminary analysis 

A preliminary analysis included the pooling of sexes and the pooling of specimens from 
all mainland Greek localities. Sexual dimorphism was tested in the largest sample from 
Peloponnisos (34 males and 49 females) and was not significant (MANOVA, p > 0.05). 

As the number of specimens collected from different localities of the mainland is 
rather low, they were pooled and tested by the inspection of the principal component 
plots (the first two, accounting for 53% of the total variability, are shown in Fig. 3). No 
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pattern of variation was found among the specimens from these localities. Similar results 
were obtained for other component plots (not shown). Thus, the pooled specimens from 
the Greek mainland will be referred to as the Continental Greek population. (CG in 
Tab. 1). 
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Fig. 3. First and second principal component (PC) plot of the continental Greek population (CG). Bor¬ 
ders of the largest subpopulations (Macedonia and Thessalia) are also plotted for comparative pur¬ 
poses. Abbreviations as in figure 1. 


Canonical analysis (CA) 

CA of Greek populations 

From the eleven Greek populations studied (Tab. 1) the first five were considered as main 
populations because their samples comprise a sufficient number of individuals and repre¬ 
sent the known subspecies of Erinaceus concolor distributed in Greece (except for those 
from the island of Lesvos). The six further populations that come from differenf islands of 
Greece all contain a small number of specimens and one cannot be certain to which sub¬ 
species they belong. The main Greek populations were subjected to CA. MANOVA was 
significant (Wilk’s A = 0.0352, F = 8.44, p < 0.0001) as well as all pairwise contrasts be¬ 
tween them (p < 0.0001). Canonical scores are plotted in figure 4. It is clear that the first 
canonical variate (CV) discriminates the insular from the mainland populations while the 
second one discriminates the Cretan population {E. c. nesiotes) from that of Rhodos 
{E. c. rhodius). The population from Lesvos of an uncertain taxonomic status, lies inbe- 
tween of all (Fig. 4). From the canonical coefficients in table 2 it can be concluded that 
CBL, ZB, BL, PL, and RA (1st CV) are the measurements playing an important role in 
the discrimination of mainland from insular populations, while UT, LT, CBL, GSL and 
PB are those that mostly discriminate Rhodos and Cretan populations (2nd CV). 
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Fig. 4. First and second canonical variate (CV) plot of the five main Greek populations (nos. 1-5). Ab¬ 
breviations as in table 1. 


Table 2. Standardized canonical coefficients for canonical analysis of the five main Greek populations 
(1-5). First CV accounts for 65.6% and second for 17.3% of the total between group variability. 



IstCV 

2nd CV 

3rd CV. 

4th CV. 

HB 

-0.046 

-0.153 

-0.306 

0.352 

TL 

0.145 

0.440 

0.321 

0.034 

EL 

0.022 

-0.112 

0.102 

-0.192 

HFL 

-0.082 

0.063 

-0.314 

-0.261 

GSL 

0.223 

0.540 

-0.664 

-0.050 

CBL 

-0.877 

-0.630 

0.352 

0.615 

BL 

0.554 

-0.344 

-0.056 

0.180 

ZB 

0.617 

-0.161 

0.470 

-0.349 

MB 

-0.388 

0.378 

0.618 

-0.054 

IC 

-0.246 

-0.096 

-0.139 

-0.215 

NL 

0.279 

-0.019 

0.334 

-0.160 

RB 

0.073 

0.206 

0.778 

0.017 

PB 

0.039 

-0.504 

-0.333 

0.136 

PW 

-0.018 

-0.025 

0.347 

0.075 

UT 

-0.378 

1.243 

-0.348 

-0.113 

ML 

-0.028 

0.179 

-0.818 

-1.221 

LT 

0.069 

-0.533 

-0.177 

1.144 

HS 

0.025 

-0.042 

-0.306 

-0.468 

PL 

-0.500 

-0.397 

0.529 

0.455 

HRM 

-0.165 

-0.442 

-0.119 

-0.065 

RA 

-0.459 

-0.073 

-0.276 

0.111 
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The insular populations (no. 6-11 in Tab. 1) were subsequently classified into the main 
classes and the results are shown in table 3. It is clear that the population from Kyclades 
(KY as well as those from the Ionian islands, i. e. Kephalonia (KP) and Kerkyra (KE), are 
all classified into the Cretan (CR) population. The population from Samos (SA) is mostly 
classified in the population from Rhodos (RH). The populations from Chios (CH) and 
Kos (KO) have uncertain positions. Note that the specimens from Kos are karyologically 
closely related to Rhodos population (Mandahl 1978). 


Table 3. Subsequent classification of the small Greek insular populations (6-11) into the main cate¬ 
gories (1-5). Abbreviations as in table 1. This classification was based on the Mahalanobis distances 
from each specimen to the centroids of the main populations. 



CR 

PE 

CR 

RH 

LE 

KY 

_ 

_ 

4 

_ 

_ 

KE 

- 

- 

5 

- 

- 

KP 

- 

- 

3 

- 

- 

KO 

2 

- 

2 

- 

1 

SA 

- 

- 

- 

5 

2 

CH 

1 

- 

2 

2 

1 


CA of the Bulgarian populations 

Canonical variate plot of the three Bulgarian subspecies (Fig. 5) described by Markov 
(1957) reveals that E. c. roumanicus and E. c. drozdovskii are clearly different subspecies 
(MANOVA, p < 0.0001), while E. c. bolkayi is an intermediate form significantly different 
from E. c. roumanicus (p = 0.005) but not from E. c. drozdovskii (p = 0.022). 


8 



istcv 


Fig. 5. First and second canonical variate (CV) plot of the three Bulgarian populations (nos. 12-14). 

Abbreviations as in table 1. 
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CA of Greek and Bulgarian populations 

Using the 15 measurements (common to all Greek and Bulgarian material, see Tab. 1) an 
8-group CA including both main Greek and Bulgarian populations was made. The 
UPGMA tree (Fig. 6) based on Mahalanobis distances revealed that the Greek mainland 
populations (CG and PE) are different from all the Bulgarian populations (MANOVA, 
p < 0.001). Insular Greek populations clearly form another cluster, with the exception of 
the island of Lesvos which is closer to the continental Greek populations. 


32 24 16 8 0 

I-1-1-1_I 

__ CG 

_ _ PE 

- LE 

- W 

_ _ Bb 

_ Br 

__ CR 

- RH 

Fig. 6. UPGMA phenogram of the five main Greek (nos. 1-5) and the three Bulgarian (nos. 12-14) po¬ 
pulations, based on the Mahalanobis distances between centroids. Abbreviations as in table 1. 


Comparison of neighbouring populations from southern Balkan Peninsula and Asia Minor 

The Yugoslavian population of hedgehogs from Vis island, described as E. c. nesiotes 
(Dulic and Tvrtkovic 1979), was compared to the Greek insular populations from the is¬ 
lands of Crete, Rhodos and Lesvos. The Mahalanobis distances calculated confirm that 
the Yugoslavian population is closer to the Cretan form than to the other Greek insular 
populations, but remains significantly different from it (p < 0.001). 

The specimens from the Ionian islands, Kerkyra and Kephalonia, were pooled in one 
population and compared to Yugoslavian insular population from Vis, because these is¬ 
lands are geographically closer to Vis than any other Greek island. We found a significant 
difference between them (p < 0.001). 

The mainland Yugoslavian populations, representing E. c. bolkayi (Dulic and Tvrtko¬ 
vic 1979) and E. c. drozdovskii (Martino 1933; Petrov 1940), due to the small size of 
samples cannot be analysed with confidence. The Mahalanobis distances between cen¬ 
troids indicated that: a) Yugoslavian bolkayi is closer to Bulgarian drozdovskii than to its 
Bulgarian counterpart, b) Yugoslavian drozdovskii is closer to its Bulgarian counterpart 
than to any other population, c) Yugoslavian bolkayi is very close to the Greek mainland 
populations from continental Greece and Peloponnisos d) Yugoslavian drozdovskii and 
Yugoslavian bolkayi are very closely related. 

Finally, the specimens from Asia Minor, described as E. c. transcaucasicus by Wett- 
STEiN (1966), were compared to all other populations. With only two measurements avail¬ 
able (CBL and ZB), they were found to be different from all populations, except those 
from Lesvos (MANOVA, p = 0.20). 
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Indices 

The basic statistics for Mi and Li in all populations are given in table 4. Length index is 
generally constant among the studied populations of Erinaceus concolor, except those 
from the islands of Rhodos and Samos which differ significantly from the others (ANO- 
VA, p < 0.001). 

Maxillary index varies amongst the five main Greek populations and the insular Yugo¬ 
slavian population (no. 1-5 and no. 17 in Tab. 1) while the ANOVA (Tukey’s test), groups 
together the populations from continental Greece, Peloponnisos, Crete and Vis island, on 
one hand, and the populations from Rhodos and Lesvos, on the other. These results are 
in accordance with our karyological data which differentiate the specimens of the Greek 
mainland and Crete from those of the islands of Rhodos and Lesvos (Giagia and On- 
DRIAS 1980). 

Table 4. Basic statistics (means and standard deviations = sd) of the two indices in the 18 populations 

studied. 


Population 

Length index 

Maxil. index 


mean 

sd 

mean 

sd 

CG 

1.66 

0.05 

1.07 

0.14 

PE 

1.69 

0.05 

1.14 

0.11 

CR 

1.67 

0.05 

1.15 

0.10 

RH 

1.56 

0.05 

0.77 

0.09 

LE 

1.64 

0.04 

0.92 

0.05 

KY 

1.69 

0.05 

1.10 

0.18 

KE 

1.64 

0.04 

0.94 

0.04 

KP 

1.63 

0.02 

0.88 

0.07 

KO 

1.63 

0.02 

0.95 

0.16 

SA 

1.58 

0.04 

0.88 

0.05 

CH 

1.61 

0.06 

0.74 

0.13 

Bd 

1.68 

0.09 



Bb 

1.66 

0.05 



Br 

1.66 

0.07 



Yd 

1.66 

0.05 

1.21 

0.21 

Yb 

1.62 

0.07 



Yn 

1.62 

0.04 

1.18 

0.11 

Tt 

1.61 

0.05 




Discussion 

The results of the above morphometric analysis of the Greek populations support the pre¬ 
viously expressed view that one large-sized subspecies of eastern hedgehog is distributed 
throughout mainland Greece (Giagia and Ondrias 1995). 

The analysis of Bulgarian hedgehogs, based on data given by Markov (1957), showed 
that two distinct subspecies occur throughout the country; the small-sized E. c. roumani- 
cus and one large-sized subspecies including bolkayi and drozdovskii, which were pre¬ 
viously described as being different. 

The comparison of the neighbouring populations from Yugoslavia, Bulgaria and 
Greece showed that previously there was confusion concerning the name of the large¬ 
sized hedgehog described by Martino (1930, 1933) and reported by Holz (1978) and 
Kratochvil (1980) as occurring in the southern Balkan peninsula (Malec and Storch 
1963; Petrov 1940; Wettstein 1942). 





138 


Eva B. Giagia-Athanasopoulou and G. Markakis 


Although the populations from the mainland of Greece and Bulgaria are quite clearly 
different, they do overlap to some extent and are characterized by morphometric continu¬ 
ity. It is suggested that one large-sized subspecies E. c. bolkayi is distributed throughout 
the southern Balkan peninsula. 

The CA of the Greek insular populations shovv^ed that iwo small-sized subspecies are 
distributed throughout Greek islands. E. c. nesiotes (Bate, 1905) from Crete appears to oc¬ 
cur also on other islands of the Aegean sea such as Syros, Tinos, Naxos in the Kyclades, 
as well as on the Ionian sea islands, Kerkyra and Kephalonia. According to Dulic and 
Tvrtkovic (1979) the distribution of E. c. nesiotes extends up the Adriatic sea to the is¬ 
land of Vis. The analysis showed that the Yugoslavian insular population is rather differ¬ 
ent and must be reexamined. The subspecies E. c. rhodius (Festa, 1929) from the island of 
Rhodos is also distributed on Samos and probably Chios, as the karyological data confirm 
(Giagia and Ondrias 1980). The hedgehogs from the island of Kos although karyologi- 
cally identical to those from Rhodos (Mandahl 1978) hold an uncertain position. De 
Beaux (1929) described some hedgehogs from the island of Kos as belonging to 
E. c. transcaucasicus which is found in Asia Minor. 

Finally, the position of hedgehogs from Lesvos island in CA and their similarity to the 
specimens from Asia Minor, support the view that E. c. transcaucasicus is probably dis¬ 
tributed on Lesvos, although the set of measurements is not sufficient to draw definite 
conclusions. This may be supported by the fact that many species distributed in Asia Min¬ 
or are also found amongst Lesvos fauna (Giagia et al. 1982; Ondrias 1966). An extensive 
study of hedgehogs from Asia Minor as well as the study of more specimens from the 
East Aegean sea islands might clarify the relationships among the insular populations 
(£■. c. rhodius) and those from Asia Minor (E. c. transcaucasicus). 
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Zusammenfassung 

Multivariate Analyse morphometrischer Merkmale beim Ostigel Erinaceus concolor in Griechenland 

und benachbarten Gebieten 

Von der siidlichen Balkanhalbinsel sind eine Reihe von Unterarten des Ostigels Erinaceus concolor be- 
schrieben worden. Trotz umfangreicher morphologischer und karyologischer Beschreibungen ist ihr ta- 
xonomischer Status immer noch unklar. In der vorliegenden Studie wird versucht, den Status der 
verschiedenen griechischen Festland- und Inselpopulationen zu identifizieren und ihre Beziehungen zu 
Subspezies in Nachbarregionen zu klaren. Dazu wurde eine multivariate Analyse von 21 somatischen 
und kranialen Merkmalen an 191 griechischen Exemplaren sowie anhand veroffentlichter Daten (76 
aus Bulgarien, 21 aus Jugoslawien und 9 von Kleinasien) durchgefiihrt. 

Die Analyse zeigt, daB auf dem gesamten Festland der siidlichen Balkanhalbinsel die groBwiichsige 
Subspezies E. c. bolkayi verbreitet ist. Zwei Insel-Unterarten mit geringerer KdrpergroBe kommen auf 
griechischen Inseln vor: E. c. nesiotes auf Kreta, den Kykladen und den lonischen Inseln, E. c. rhodius 
auf den ostagaischen Inseln Rhodos, Samos und Chios. Die Subspezies E. c. transaucasicus kommt auf 
den Inseln Lesbos und Kos vor. 
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